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This paper is devoted to a detailed investigation of the perturbed pendulum-like motions of a 
heavy rigid body about a fixed point. Canonical variables that allow one to simplify the analysis of 
homoclinic and heteroclinic orbits are introduced. Characteristic properties of perturbed pendulum- 
like motions of the body in inertial space are studied. A qualitative description of asymptotics of 
pendulum-like motions in a neighbourhood of split separatrices is given. 
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(N ■ I. INTRODUCTION 

t-H ■ 

Differential equations of motion, referred to the principal axes inertia of the body at its fixed point, have the 
\ following form 

— = M x u + /i(7 x r), ^ = 7 xu, (1) 

where uj is the angular velocity, M = Au> is the angular momentum of the rotated body, A = diag(Ai, A2, A3) is the 
tensor of inertia, 7 is the unit vector directed vertically upward, r is the unit vector to the body's center of gravity, 
(i is a product of the weight and the distance from a fixed point to the gravity center. 

Let us start the investigation of a body motion, whose mass distribution is restricted by the following conditions 

A, > A 3 > A 2 , n > 0, r 2 > 0, r 3 = 0, ft < 1. (2) 
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If the conditions ([2]) have been fulfilled, and the following equalities wi = ui 2 = 73 = are used as the initial values 
of the variables, so the body is rotated about the fixed horizontal axis as a physical pendulum: 



Ml = M 2 = 73 = 0, M 3 = ± y / 2A 3 (ho +70 dnr, 

r\ [2 sn 2 r — 1] =p 2 T2 snr cnr, 72 = r 2 [2 sn 2 r — 1] ± 2 n snr cnr, 



71 



(3) 



(t-t )VOo + M)/(2A 3 ), k = j2n/(ho + n), 



where the energy constant ft-o > M is defined by the initial conditions. 

For the further investigation of a rigid body motions in a neighborhood of the integrable Euler case the canonical 
Andoyer-Deprit variables are successfully used pH^. We will use the following notations: G is the modulus of the 
angular momentum M, L is a projection of the angular momentum M on a moving axis Oz, H is a projection of 
the vector M on an upward-directed vertical; the variables canonically conjugate to G,L,H are the angles g,l,h 
in range [0, 2ir). The dependency between phase variables and Andoyer-Deprit 's ones are expressed in the following 
formulas: 



Mi = VG 2 - L 2 sin Z, Mi = \/ G 2 - L 2 cosl, M 3 = L, 



A) 



71 = (s'mqcosp + s'mpcosqcosg) sinl + cos g sing cos 

72 = (sing cos p + s'mpcosqcosg) cosl — cos q sin g sin I , 
73= sinqsinp — cospcosgcosg, 

where sinp = L/G, sing = H/G. The inverse transformation is given by 



(5) 



G = |M|, g = arcsin 



(M271 - Mi 72) 



IMI 



y/M? + M| |Mx 7 |_ 
L = M a , I = arctan(Mi/M 2 ), H = M ■ 7. 
We introduce new parameters ao,/3o,7o that characterize the mass distribution: 



(6) 



ao = 2 arctan 



To 



2 arctan 



a/Ai(A 3 - A 2 ) 
y/A 2 (A 1 - A 3 ) 



WA 3 -A 2 



f3 — 2 arctan 



e [o.tt] 



(7) 



.-/Si - A3. 

Immediately from ([7]) the following relations can be found 

,2 "o „o 



where 7r > ao > 70 > 0. 



sm 



Ai = A 



.2 7o 



A2 = As 



2 ^° 

3 y 



'"1 



cos ■ 



Pa 



'2 



sm ■ 



If ao = /3q, then the rigid body mass distribution is subjected to the Hess conditions (the center of gravity belongs 
to the perpendicular, drawn from a fixed point to the circular section of gyration ellipsoid). If 70 = /3q, then body 
mass distribution is subjected to the Grioli conditions (the center of gravity belongs to the perpendicular, drawn from 
a fixed point to the circular section of inertia ellipsoid). Moreover, from the triangle inequality A 2 + A3 > Ai we 
determine an additional restriction for the parameters: 



cos ■ 



a 



With regard to (J5]) the Hamiltonian of the mechanical system will have the following form 



L 2 )(A^ sin 2 I + AT 1 cos 2 I) + A^L 2 ] 



+H [singcospsin (I + /3q/2) + cosgsinpsin (I + (3q/2) cos <?+ 
+ cosqcos (I + Pq/ 2) smg] . 



(8) 



We changed the sequence of the principal inertia axes in moving basis Oxyz, therefore thepresent canonical variables 
G,L,H, g,l,h are distinct from the standard variables, used in the research works [iMa]- For rotations of a rigid 
body about its intermediate inertia axis, the canonical variables ([6]) are degenerated. Later, it will be found, that 
such variables simplify the study of the perturbed pendulum motions. Under conditions @ the rotations ([3]), as it 
is known, are unstable. Due to this, the investigation of qualitative properties of the perturbed body motions in the 
fixed space is of our special interest and also determines the main objective of this paper. 
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II. NEW CANONICAL VARIABLES 



Using the generating function 

S 2 = (21 + 0o)Ji + (g~l- Aj/2) Ji + hJ 3 

we shall find the relations 

L = -7TT- = 2Ji — J-2, G = — — = J2, H = -7-— = J3, 

dl dg dh ^ 

aJi oJ 2 oJ 3 

which define the canonical transformation to the variables (Jj, 0j), i = 1,3. So from ([9]) the following expressions are 
obtained 

! = ^, 5=^1+^2, ^ = 3l Ji = ^±A J 2 = G, J 3 =if. 



Let us introduce the dimensionless variables J[ — Ji^J(A\ — A^f/jAxA^)-, i = 1,3, and the notation for time t' = 
ty/(Ai — A 2 )/(AiA 2 ). In order to reduce the record, the primes of dimensionless variables should be omitted. In the 
new variables the Hamiltonian ([SJ is transformed to the following form 

H =Ji{J% - J\) [cos (6>i - p Q ) - cosa ] + i xjf + 

/ j2 _ j2 

+/i [Ji sin (6»i + 2 ) + (J 2 - Ji) sin0 2 ] — 1 + (10) 
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+ 2 -" 72 sm y ' 

where 

A\A 2 _ sin 2 «o 



> 0. 



A 3 (Ai~A 2 ) 2 (cos 70 - cos ao 
At J 3 = level the Hamiltonian (ITU|) will be expressed as follows 

H = ^ Jf + Ji( J 2 - Ji) [cos (6»i - /3 ) - cos a ] + 

t (11) 
+4- [Ji sin (0i + 2 ) + (J 2 - Ji) sin 2 ] , 

We will suppose that angle 8 2 is changed in range [0, 2n), and the alteration of variables Ji,2 is restricted by the 
following inequalities J 2 > 0, J 2 > J% > 0. At J 3 = level, the differential equations of rigid body motion are the 
Hamilton equations 

4 _m f _ _m 

l ~ dJC l ~ 90,' ' ' 

we obtain them in the explicit form: 

0'i = (J 2 - 2Ji) [cos (0! - fa) - cosa ] + [sin (6 1 + 2 ) - sin0 2 ] , 

A = Ji(J 2 - Ji) sin (0x - p ) - ^ cos (0! + 2 ), 

Ji ( 12 ) 
2 = x,J 2 + Ji [cos (0i - /3 ) - cosa ] - fJ>~j2 [sin (0i + 02) - sin0 2 ] , 

J2 

J2 = — t [Ji cos (01 + 62) + (J 2 - Jl) COS0 2 ] . 
J2 

If we know the system (|12p solution, so the cyclic coordinate 3 can be easily defined by the quadrature from the 
equation 



% = 2fi -3 sin-. (13) 
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III. PARTICULAR SOLUTIONS OF THE HAMILTONIAN SYSTEM ([12]) 

For the small values of the parameter /i, the angular momentum modulus |M| = J 2 {t) remains in a small neigh- 
borhood of its initial value J 2 = J 2 (to). Let us consider the exact solutions of the system (JT2J) , which correspond to 
the energy levels, close to h = xJ 2 (t )/2 > fj,. 

A. Unperturbed system solutions 

As fi = we have the unperturbed system of differential equations 

01 = (J 2 - 2 Ji) [cos (0i - P ) - cosa ] , 

2 = xJ 2 + J\ [cos (0i - p ) - cosa ] , (14) 
ii = Ji{-h - Ji) sin (0i - p ), Ji = 0. 

This system is completely integrable: it describes the rigid body motion in the integrable Euler case. The general 
solution of the system (| 14[) may be expressed in terms of theta functions of time. The unperturbed separatrices, used 
below, correspond to doubly- asymptotic solutions of the system (fl4l : 



1) 0i = Po - «o, h = 4/ (> (t -* o)+w + 1) , 

2 = xJ 2 o (i-i o ) + 0°, J 2 = J 2 °; 

2) 0i = a + po, Ji = J 2 / ( e -«(*-*o)+* + l) , 
2 = x4(t-t )+el J 2 = J 2 °; 

3) 0i = /3 - «o + 27T, Ji = J 2 / ( e «( t -*«)+- + l) 

2 = ^j 2 °(<-< o ) + 02, ^2 = J 2 °; 



(15) 



(16) 



(17) 



4) 0i = a + po + 2tt, Ji = J 2 / ( e -»(*-*o)+» + 1) 

e 2 = Kj°{t-t ) + e°,j 2 = j 2 °, 



(18) 

9 2 = ^{t-t Q )+9 v 2 ,J 2 = J 2 U , 
where 0^, Jj are the initial values of the variables at t = to, the constant parameters are denoted through 



(jo _ JO) 

u = J 2 sina o >0, u = In — 2 „ 1 € (—00, 00). 

•h 

In the limiting cases from the relations (fT5|) — (fT5|l the following equalities can be found: 

1) 9i=p Q - ao, Ji = 0; 1*) 0i = ft, - a , J x = J 2 ; 

2) 6i=p + ao, Ji - 0; 2*) 0i = p + a , J x = J 2 ; 

3) 6 1 =Po-a + 2tt, Ji = 0; 3*) 0i = Po - a + 2tt, Ji = J 2 ; 

4) 9 1 =p +a + 2tt, Ji = 0; 4*) 0i = /3 + «o + 2tt, J x = J 2 . 

At Fig. [T]the phase orbits of unperturbed system of equations (fT4")) at the sphere |M| 2 = const and at the plane 
K 2 (0i, Ji/J 2 ) for the values A\ = 2, A 2 = 1, A 3 = 1.5 can be reviewed. The separatrix, composed of two intersecting 
circles, Fig. [TJ a, is transformed in the vertical segments, Fig. [TJ b. 



Figure 1: Phase orbits of the Euler case: a) at the sphere |M| 2 = const, b) at the plane R 2 (#i, J1/J2). 



B. Physical pendulum 

Pendulum-like rotations of a rigid body are expressed explicitly by Jacobi elliptic functions of time: 
1) J 1 =0, J 2 = ^J2(h + dnr, sin6> 2 = 2sn 2 r - 1, 



2) Ji = J 2 , J 2 = y/2{h + n)/x dnr, sin [6 1 + 6 2 ) = 2snV - 1, 



(19) 



where r = (t — t Q )y/ (h + Jl)x/2, k = \j2[ij (h + [i). 

From the relations (j4]),((9]) it follows, that the angle 6\ hasn't been defined for the equalities J\ — and J\ = J 2 . 
This degeneracy is not essential, as far as the body's position and velocity, relevant to the mentioned equalities, are 
uniquely defined by the smaller numbers of the phase variables. 

If Ji = o(/Lt), then by using the integral (ITT]) we exclude the value of J 2 from the 1st and 3rd equations (|T2|) . Due 
to the variables transformation y = ctg(#i/2), r = 62 , we can obtain Riccati differential equation 

y' + a y 2 + ai {T)y + a 2 {r) = 0, (20) 

1 r n \ , \ 1 • n flCOST 

a Q = — (cos/? - cosa ), Oi(rJ = — sm/3 1 



IH >f 2{h — /ismr) 

1 / « \ ^sinr 

a 2 (r) = — — (cosa + cos/? ) - 



2>tf 2{h — /isinr) 

If J 2 — Ji = o(/i), then the similar Riccati equation (|2T)|) will be obtained after transformation y = ctg(#i/2), 
t = 61+62- The dependence y(r) allows investigating 61(62) in the arbitrarily small neighborhood of pendulum 
motions. Also we can note that an analytical study of the variational equations for pendulum motions of a rigid body, 
and the simplest cases analysis of its integrability are presented in the paper [|| . 

Small perturbations of the exact solutions, considered in the cases a),b), lead to the formation of asymptotic 
solutions at the energy levels, close to ho = xjf (to)/2 > /-*• Let us introduce formulas describing the asymptotic 
solutions of the equations (1121) to the first-order approximation. 



C. The solutions, close to pendulum type 

We get the approximation formulas for small values \x. Let us put 8*,(t) = xJ 2 (t — to) + 6*2- The periodic solutions 
of the system (fT2")) are the following: 

Ti : 61 — /3 Q + ea + (J,(du s'md*. + d 12 cos #3), Ji=o(//), 

62 = 02 + cos9 * 2, j2 = J ° ~ ^ sinf?2 * ; 

T 2 : 6x = (3 + ea + fi{d 2 i sin 6* 2 + d 22 cos 62), (21) 
Ji = J 2 = J 2 ° - -^0 sin (A> + £a o + ^2). 



62 = 9% - M(^2i sin 6*2 + d 22 cos 05) H /i n cos 05, 

x(J 2 u ) 2 
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where e = ±1 and constant coefficients dij are expressed as 

<iii 



121 



x sin 


(0n - 


+- £ao) + £ sin q;o [cos (0o - 


f- £an) — 11 






(J 2 °) 2 (x 2 + sin 2 ao) 




x[l - 

^ L 


- COS 


(A) H - £<^o)] H~ ^ sin o.q sin 


f/3n 4- £(y.n) 






(J 2 °) 2 (x 2 +sin 2 a ) 




>rsin 


Wo- 


f eao ) — £ sin a [cos (0q - 


\- £a ) - 1] 






(J 2 °) 2 (x 2 +sin 2 a ) 




x[l - 


- COS 


(/?o + eao)] — £ sin ao sin 


(0o + £a ) 


(J 2 °) 2 (x 2 + sin 2 a ) 



C?22 = 



The system (|12[) solutions, relevant to the asymptotics of pendulum-like motions, come arbitrarily close to the limit 
cycles Ti, Ta, furthermore the variable Q\ is in the range \0o + £ao — Ai, 0o + £ao + Ai], where 



After calculations we obtain 



Ai = 



Ai = nyd 2 n +d 2 12 = n\Jd 2 2l + d 2 2 . 

/i sinao sin ||/3o + £ao| 
y^sin 2 a + 16 sin 2 \ (a + j ) sin 2 ±(a - 70) /io 

n Ai ^^2(^3 - A 2 ) + £r 2 A 2 y /A 1 (A 1 - A3) 



(22) 



^4 3 (yli - A2)(4AiA3 + 4A 2 A 3 - ZA^ - AA 2 ) h 



The same arguments may be used for investigation of the limit cycles, received from (|21[) , with the replacement of 
00 + sa by O + + 27r. 



D. Perturbed separatrices 

As an example, let us consider one of possible variants. For this purpose, the following expressions 

e l = a + 00 + se u Ji = j*{t) + 5Ji, e 2 = e* 2 {t) + se 2 , j 2 = j 2 ° + sj 2 

should be assumed in the neighborhood of the separatrix (TTBl . where the functions J*(t), 2 (t) are defined by the 
formulas (|16|) . The linearized equations of the perturbed motion take the following form: 

59x = (2 J* - J 2 °) sina + -4 [sin (a + 0o + 0* 2 ) - sin^] , 

•h 

SJi = (J 2 ° - 2J X *) sina 5Ji + J{{Jt ~ J{) cosa 56 x + 

J* 

+ J* sina 5J 2 - fi cos (a + 0o + 9 2 ), 



9 2 = x SJ 2 — J{ sin ao S61 + /1 



J* 



sinf 



i(a + y5 + ^)] 



Ji 



8J 2 = —fi cos 9 2 + fi — q- [cos 6* 2 — cos (a + /3 + 6*2)] 



4 

The equations (f2"3")l admit the first integral, depending on (56*1 , SJ 2 



J* 

J = J*{J* - J°) sina 86\ + >cJ^5J 2 + fi sm6 2 + fj, -4 [sin (a + /? + #2) - sin 6 2 ] 

Jo 



(23) 
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Differential equations (1231) are integrated in quadratures. In particular, near the unperturbed separatrix the small 
deviation of 9\ from the initial value 9\ — olq + /3o we obtain from the first equation of i 



Jt{j° 2 -j{) 



I J*i{J°2 - J*i) cos (9* 2 + ~(c*o + A))) dt = 

Jt 1 



( e Mt-t ) _|_ e v\2 ^ ( e «(*-*o) _|_ e v\2 

U-^ + iy ^ (<o) + 2 " Sin 2 (ao + A > J°e^V X 



t— to gUT 



(e UT + e v ) 



— cos (xJ 2 °r +-(a + M + 0°) dr. 



(24) 



IV. MELNIKOV'S INTEGRAL 



According to the notations [H, we write the Hamiltonian ([Tl~T) as sum H = Hq + f/H\. The distance between 
the stable and unstable separatrices of the system (jT2"j) may be investigated analytically by means of the Melnikov's 
integral Q 



Tm(6) 



T~Lq, -^Jq \ (J*; #|'. ■!-;>■ W o + W if / /. 



^0) "^Jo ^ (J*'®!-* ^2^2 



dH dUi dU Q dUi 



80! dJi dJi 80i 



evaluated along an unperturbed orbit (|16p . which connects hyperbolic periodic orbits. The functions J*(t), 9®, 
J 2 , #2 W correspond to the doubly asymptotic solution (|16[) . Owing to the transformations we find 



Im(0) = <to 



{e 1 -' 



■ cos (xJ°t + -(a + /3 ) + 9)dr, 



in which <to = 2>f~ 1 sinao sin i(a + A)) is nonzero constant. The integral Im{9) is a divergent improper integral, 
so we will be interested only in the evaluation of its principal value (the fast oscillating part is excluded). Using the 
linear substitution r = (2ttx + v)/u, we take into account the notations 



- >0, 5= _(ao +£„) + - , 

sm ao 2 sin ao 



The principal value Im{9) is reduced to Legendre improper integral, which is easily calculated: 



T M (9) = ^ C os(5 + 
u 



l-4u 



sin (2vx) 



dx 



e 2TTX _|_ I 

(t q v cos (6 + 9) _ 2ir sin i(a + /?o) cos (5 + i 



(25) 



sinh v 



7° 

■ J 2 



sm ao 



sinh ^ 



Therefore, for any values of parameters, restricted by the conditions ©, the integral Im(9), considered as a function 
of the real argument 9, has simple zeros only at the points 9 = ±kir/2 — 5, k — 1,3, ... . These values define two 
heteroclinic solutions of Hamiltonian system (|12l) . asymptotically tending (if t — > ±oo) to two different periodic 
solutions Ti,T2. The transversal intersection of the perturbed separatrices means that for the rigid body, that 
satisfies the conditions ([2]), chaotic motions always exist (at least, near the separatrices), if only the parameter /j ^ 
is sufficiently small in comparison to the energy constant ho. 

Following the methodology developed by Henri Poincare, it is possible to prove, that for fixed values of the pa- 
rameters ([2]) and initial conditions, which correspond to J3 = 0, ho ^ \x > 0, Hamiltonian equations (IT21) admit a 
countable set of heteroclinic solutions. These solutions describe the doubly asymptotic pendulum-like motions of a 
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1* 


f « * A 


* -i* 

















Figure 2: Heteroclinic cycle at the plane ffi 2 (#i, J1/J2). 

rigid body. A heteroclinic cycle of the system (fT2|) is shown schematically at Fig. [2j it is composed of hyperbolic 
periodic orbits and heteroclinic orbits, which lie in intersection of perturbed separatrices. 

An improper integral Im (0) is a special case of an integral evaluated by Ziglin 3] . In order to check this statement 
it is possible to use the simplified formulas, obtained by Dovbysh in [3, p. 367]. The calculations, performed in this 
paper, confirm the well-known result of Kozlov 0, p. 104] about splitting of separatrices of asymmetric body, rotating 
about a fixed point in a weak field of gravity. 

V. ASYMPTOTICS OF PENDULUM-LIKE MOTIONS 

The dynamical system ([l} under the restrictions ([2]) has two hyperbolic periodic solutions T\ and T2, which are 
represented by two closed nonintersecting curves in the space R 6 (M, 7). Two asymptotic surfaces .Sf", cross each 
of these curves. If fi = these surfaces coincide pairwise, namely S± = 5^, = . The Hess solution belongs to 
the double separatrix = S*^ which, as known from |2|, § 4.6], is not split under perturbation (if fj, ^ 0). For the 
asymptotic motions of a rigid body in the case ho > /x the following properties are fulfilled: 

• The invariant tori, which could isolate one of the following separatrices , S2 , S± , , are not existed; 

• For fixed values of parameters of the dynamical system ([!} at any energy level ho > \i there is a countable set 
of heteroclinic solutions from (1 5£j [S^ f~l S^) between Lyapunov periodic orbits T\, T 2 ; 

• Any two trajectories on the separatrix (S2 ) are separated by heteroclinic orbit, nearby trajectories diverge 
exponentially with time; 

• Various types of homoclinic and heteroclinic orbits are characterized by a finite binary sequence z(9), coding 
the route along the unperturbed separatrices (I15|) -([T8 |) . 

Projections of phase trajectories of the perturbed system (|12[) in R 2 (9i, J1/J2) are shown at Fig.[3J The calculations 
were carried out for the following values of the parameters: 

/i = 0.01, fa = 70, el = a + A), J 2 ° - 0§ € (0, 2w); 

the values of cvq, 70 can be found in Table HJ At the lower part of Fig. [3] the surfaces S± , formed by the solutions of 





a 


b 


c 


d 


e 


/ 


g 




2.91891 


2.23654 


1.23096 


1.61146 


1.97724 


0.79020 


0.56207 


70 


2.46192 


1.95519 


0.92730 


1.34948 


1.77215 


0.57351 


0.40272 



Table I: The values ao,7o for Fig. [3] 

the system (|12p . are shown, which in case t — > — 00 are asymptotically approached to the periodic solution T\. In the 
upper part of Fig.|3]the surfaces S£, formed by the solutions of the system p^|) . are shown, which in case t — > +00 are 
asymptotically approached to the periodic solution T 2 . The set So = n 5^) U (5^ n S^) consists of heteroclinic 
solutions of the system p^|) . Homoclinic solutions of the system (TT2"j) belong to the set S* = (S^ n S^) U (S2 H )■ 
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a b c d e f g 

Figure 3: Asymptotic surfaces S± ,S£ ■ 




Figure 4: The sequence z(9k) for the orbits on S 1 . 

At Fig. |4] one can see the possible routes of trajectories on 5-f in the vicinity of the unperturbed separatrices (fT5|l - 
(HHJ). Using the points 9k — k/10, where k = 0,62, we define a partition of the interval [0, 2tt), in which Melnikov's 
function (|25p argument 9 is changed. Let us consider the set of solutions for the system (|12l) . subject to the initial 
conditions 9\ = a + /3 , J® ~ 0, Q\ = Of., J® = const, in a small neighborhood of periodic solution T\. Near 
the solution T 2 , according to (|24|). the phase orbit leaves the neighborhood of the separatrix (|16p and, furthermore, 
moves along the perimeter of the left or right rectangle at Fig. [2j Let Zi(9k) = 0, if the trajectory is displaced on 
the left from the incoming separatrix, and Zi{9k) = 1, if the trajectory is displaced on the right from the incoming 
separatrix, Fig. [5] i is the index of the sequence of hyperbolic periodic solutions, in the neighborhood of which the 
solution of the system (fl"2)) under consideration has been running across. The columns at Fig. [4] are composed of 
the first members of the sequences z(#fc), where black and white cells coincide with and 1 of Zi(9k),i £ Z. The 
calculations have been rendered for i = 1, 15 and the following values of the parameters: A\ —2, — 1.0, A% = 1.5, 
r = (V2/2, V2/2,0),n = 0.001, ho = 1.3342. The bottom row of the array at Fig. H corresponds to i = 1, the top 
row is that of i = 15. Black and white cells at the bottom row, Fig. SI correspond to different signs of Melnikov's 
function (|25[) . Right array of Fig. 0] is a magnification of some part of the left array: using the following points 
9k = 37r/40 + fc/100, where k — 0,31, the partition of the interval [37r/40, 77r/40] is defined. Multiple zooming of a 
partition step has no effect on the general qualitative properties of the phase orbits under investigation. 
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VI. BODY MOTION IN A FIXED BASIS 

We will study the qualitative properties of body motion in a fixed basis. With zero constant of the momentum 
integral (J 3 = 0) the angular momentum M during all time of motion, lies in a horizontal plane OrjC,. The direction 
of vector M in a plane OrjQ and its modulus |M| are characterized by variables 9 3 , J 2l therefore the time evolution 
of the angular momentum is described by the functions 63(f), ^(f)- It is well known that the pendulum motions, 
which are described by (JTSJ) , must satisfy the additional restriction 3 — const, i.e. the direction of M is preserved in 
a fixed basis, but the modulus |M| is changed. For the perturbed separatrices, described by the linearized equations 
of type (|23p . it is possible to find the increment of angle 63 with increasing time to infinity: 

(1) 2/i7r sin±(/3 -ao) fiirx sin ±(/3 - a ) 

1) 9 1 = a - p a , lA y >9 3 » » ; » — : , 

{J 2) sieqo ilq smao 



(2) 2 fin sinjffo + ao) finx sin + a ) 

2) V x = a + Pq, L\ y 'V 3 w - » : w — : , 

(0/2 J smao "0 smao 

3) 0? = a - A) + 2tt, A^ = -A«0 3 , 

4) 6? = a + /So + 2tt, A< 4 »e 3 = -A< 2 >0 3 . 

In particular, for Hess case (ao = /3o) in paper 9] the following expression has been obtained 



(26) 



A (2)( 



2/Z7T fnrAiA 



2 



(J 2 ) 2 A 3 {A X -A 2 )h 



It should be noted, that this value can be large enough, if the inertia ellipsoid is slightly different from the sphere. 
For Grioli case (70 = Po) from (|2"6")l , including (JTJ) we obtain 

AWfc a ^ Slna ° e f ^ (1 - V2) \ ^ 
4/i sin±(a +70) I h y/2 'J' 

4/i sin 2 (a - 70) 

The angle 63 is the angle of precession of angular momentum about the vertical axis. In space of rigid body 
parameters we can distinguish subsets with different angle 9 3 behavior in the neighborhood of perturbed separatrices: 

i) 7T = p > a > 0, A< 2 >0 3 = A«0 3 > > A< 3 >0 3 = A* 4 '^; 

m) tt > A > a > 0, A^0 3 > AW03 > > A( 3) 3 > A< 4 'fl 3 ; 

in) 7T > p = a > 0, A( 2) 6» 3 > A^Ss = = A( 3) 3 > A^fig; (27) 

£«) 7T > a > /S > 0, A (2) 3 > A (3) 3 > > A (1) 3 > A (4) 3 ; 

0) 7T > a > A, = 0, A (2) 3 = A (3) 3 > > A (1) 3 = A (4) 3 . 

Thus, for small values of y, the perturbed motion orbit returns to a small neighborhood of separatrices (|15p - (|18[) . in 
this case the angle 9 3 gets the increment A^0 3 in agreement with formulas (f2l)|) . At the other areas of motion near 
boundaries J\ = 0, J2 — J\ = we have A03 = o (/z). The distinctive feature of heteroclinic and homoclinic solutions 
is the existence of finite limits 

lim 3 (i) = 67, lim 3 (t) = 0+ 

In general case, doubly asymptotic solutions orbits are considered to be the heteroclinic in a fixed space, but it is 
possible to select the dynamical system parameters in such a way, that the following difference (9 3 — 97 ) ^ could 
be any value, e.g., a rational multiple of 2-7T. 

As a geometric interpretation of a rigid body motion let us consider Poinsot kinematic representation by the rolling 
motion of the bodys ellipsoid of inertia on a fixed plane in space. Let us /1 = 0. The inertia ellipsoid for the fixed 
point O is defined by the equation 

Aix 2 + A 2 y 2 + A 3 z 2 = 1. (28) 
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We denote by p value the radius-vector of this ellipsoid, directed on the instantaneous axis of rotation, and draw the 
tangent plane IT through the end of the indicated radius- vector. Then the inertia ellipsoid ([25} during the motion 
rolls without sliding on one of its tangent planes, this plane II is orthogonal to the angular momentum M and remains 
fixed in the space. The angular velocity u> of a moving body is directed along the radius- vector of a point of contact, 
and its value is proportional to p. Louis Poinsot developed a visualization to the motion of the endpoint of the 
angular velocity vector. The path traced out on the inertia ellipsoid by the angular velocity vector u> is called the 
polhode; the corresponding curve on the invariable plane is called the herpolhode. Geometric properties of polhodes 
and herpolhodes essentially depend on the values of moments of inertia and initial conditions of motion. 

The polhode is a closed curve encircling the major or minor axis of inertia ellipsoid. The boundary curve, which 
splits two families of polhodes, consists of two intersecting ellipses, passing through the intermediate axis of inertia. In 
this case |M| 2 = 2A 3 h , then the polhode is called the separatrix. Herpolhode, corresponding to polhode-separatrix, 
is a symmetric spiral in the fixed plane IT, which runs around the Q - an intersection point of the momentum M with 
the plane IT. The length of this spiral is finite, and equals to the corresponding polhode arc. 
















|d.i |\(o 


-\Aty®rFy J ; ) 2 / 







a) a(t) for B\ = 1.02 




b) (w£(t),w f (i)) for 6% = 1.02 




c) a(t) for 6% = 1.83 d) (u((t), U((t)) for 8% = 1.83 

Figure 5: Perturbed herpolhodes in a neighborhood of pendulum-like motions. 



For p ^ we write the equations of perturbed Poinsot herpolhode. For this purpose it is necessary to find the 
components of angular velocity vector w = uigii + w T) i2 + ujq iz in a fixed orthonormal basis (ii, i2, is). 

Assuming, that oj v = oleosa, u)£ = w p sina, we can use the Kharlamov's [l(| kinematic equations describing the 
fixed hodograph of an angular velocity: 



uj£ = oj ■ -y, 



H 2 -(^-7) s 



2 da 
Jp ~dt 



(7 x w) 



~dt' 



(29) 



For the pendulum motions (|3|) let us assume, that cj^ = = 0. Under a small perturbation of the pendulum rotations 
the endpoint of the angular velocity vector traces out the curve (wj (t), Wf(£)) in the fixed space, its projection 

on the plane R 2 (w^,a;^) is a spiral: a planar curve that winds around the origin of coordinates. Near the separatrices 
(IT5|) - lfT8|) the angle a(t) receives the nonzero increment, defined by formulas (|26|) . Indeed, under the restriction J3 = 
the angles #3, a are connected by the following relation 



i(t) = a{t) — arctan 



/M ■ (u x 7) 

M-UJ 



(30) 



At Fig. [5] for two values of 6\ the graphs of the function a(t) and the curves (u)^(i),ui^(t)) are presented. These 
calculations were carried out for the following values of the parameters: 

A= (2,1.5,1.8), r= (0,1,0), /u = 0.01, h = 1.1111. 

In this case the following conditions (|27|).z): A^ 1 )^ = A< 2 )0 3 = -A( 3 '6» 3 = -A^6» 3 w 0.057715 have been fulfilled. 
The motion along separatrices (IT5"1) - (|18[) . as it follows from the results of the numerical integration of the system (fTJ 
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may be described by the ordered sequence of the heteroclinic cycle tops at Fig. [2j 

6% = 1.02 : {2 - 2* - 1* - 1 - 4 - 4* - 3* - 3 - 4 - 4* - 1* - 1}; 
0l = 1.83 : {2 - 2* - 3* - 3 - 4 - 4* - 1* - 1 - 2 - 2* - 1* - 1}. 

The change in the angle a(t) depends on the route of motion along unperturbed separatrices (see Fig. [5Ja,c). At 
Fig. GUM the curves (td^ft), u>^{t)) are shown, which coincide with the different routes of motion from the top 2 to the 
top 4* of heteroclinic cycle. 
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